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ABSTRACT 


/ 


A code has. been developed for the calculation of interplanetary solar 
proton 1 1 uences at 1 A. U. for the active years 1977-1983. The fluences 
are presented as functions of mission duration t* energy threshold E, and 
confidence level Q. For a given combination of x and Q* the routine deter- 
mines whether ordinary or anomalously large events are to be considered, 
and in the latter case, the number of anomalously large events that are 
predicted by probabilistic theory for the specified mission duration. The 
code is described in detail. A listing and sample calculations are 
attached. 
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SOLPRO: A Code to Cal culate F.nergetlc Solar Proton Fluences 

Introduction 

The code was originally developed as a supplement to UNIFLUX, a unified 
orb'tal flux integration system ( Stassinopoulos and Gregory , 1975), from an 
initial solar proton model ( Stassinopoulos and KlnR , 1979) designed exclu- 
sively for the evaluation of flux levels to be encountered during the 
active phase of the next solar cycle (1977-1983) by ear th-oroi ting space- 
craft on missions involving partial magnetospheric shielding. 

This first empirical model was characteristically simple and easy to 
use. It was derived entirely from experimental satellite measurements 
made during the active years of solar cycle 2C; it included only crude 
confidence level estimates, that is, the probability of actual cycle 21 
fluxes exceeding the predicted intensities. 

A more probabilistic analysis of solar cycle 20 fluxes was later given 
by King (1979). This formal statistical treatment was subsequently used 
by King and Stassinopoulos (1975) in their computation of ratios of solar 
to t'- ipped proton fluences for circular-orbit geocentric space missions 
to be flown during the active years of the next solar cycle (1977-1983). 

In this more sophisticated model the probability of exceeding mission 
fluence levels is given as a function of mission duration and energy. 

The new SOLPRO code incorporates all the important features of the 
King (1979) analysis, which includes the distinction between "ordinary" 

(OR) and "anomalously large" (AL) events, and the probability of occur- 
rence of the latter. Be it noted, that only one anomalously large event 
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was observed during solar cycle 20: the event of August, 1972; according 
to King (1974), the proton fluences of that one event alone constituted 
about 697. (E >10 MeV) and 847. (E>60 MeV) of the entire cycle fluence. 
Consequently, whenever at least one AL event is predicted, the OR-event- 
contribution is negligible. 

The referenced solar proton models arc based on the most continuous 
set of satellite data available to date, covering the period of enhanced 
solar activity during the 20th solar cycle (1964-1975) through the 1972 
AL event, after which no significant events occurred. In both models, 
all fluxes were taken to be isotropic. As noted in Stassinop oulo s and King 
(1974), on an event-integrated basis, departures from isotropy are typically 
only a few percent or less. This applies to the interplanetary medium 
and to much of the magnetosphere, but may not apply at low altitudes 
where anisotropy should result from atmospheric loss mechanisms. 

When using the SOLPRO routine, two important items should be kept 
in mind. First, there is no assurance that the overall flux levels 
observed during the 20th solar cycle will also occur during the 21st 
cycle. However, as indicated in King (1974), a comparison of the annual 
mean sunspot numbers of past solar cycles suggests that cycle 21 will most 
likely be similar to the very ordinary cycle 20 than to the very extra- 
ordinary cycle 19. Second, there is no reliable way of predicting the 
distribution of individual solar events in time, in flux level, and in 
spectra through the 2lst cycle. That pertains to ordinary as well as to 
anomalously large events. Therefore, the results obtained from the routine, 
although the best predictions available of solar proton fluences expected 
to occur during 1977-1983, should be considered guideline approximations only. 
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In the code, the solar proton fluences are expressed as functions 
of the parameters "confidence level", Q* (in percent), "mission duration", 

T (in months), and energy, E (in MeV). Within the limits established by 
the routine the user is free to select his own values of Q and ^ for pro- 
cessing, but E is not an input parameter and the calculations are performed 
at ten fixed energy thresholds. 

For missions whose trajectories will involve significant time away 
from a heliocentric distance of 1 A.U. (astronomical unit), a helioradial 
dependence of event fluences must be allowed for, as discussed by Kin^ (1974). 

Similarly, for missions whose trajectories will involve partial 
magnetospher ic shielding, the fractional exposure factor has to be ' 
determined, as discussed by Stasslnopoulos and K ing (1974). An independent 
code performing these shielding calculations has been developed (Stas ■; inopouios , 
1975) and is available from the National Space Science Data Center, Grcenbelt, 
Mary land. 

A detailed descr iptl on of the solar proton routine Is given in the 
next section on "Method and Procedure", with an analysis of the method 
employed and a review of the organization and structure of the code. In 
section "Application and Use" the routine is briefly discussed, its opera- 
tion is described, and the arguments of its transfer vector are presented. 

Code listing and sample outputs are given in the attachments. 


^ - 1 ' — - ' — ■ ■ - - .... • . . - . ■ . . . 

Q is the degree of confidence a mission planner wishes to have that the 

actuutti/ zncou.ntc.fLzd solar proton fluence will not exceed the 
pfLzd-ictcd fluence, as given by the statistical solar proton model used. 
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Method and Procedure 


SOLPRO is a short, compact routine which evaluates omnidirectional, 
unattenuated, interplanetary energetic solar proton fluences in terms of 
three variables: Q, T, and E. The first two are input parameters, the 

last is internally programmed. As defined elsewhere in this report, Q 
denotes the degree of confidence one wishes to assign to the results, 
namely that for the specified mission duration t the calculated fluences 
at ten preset energy values ranging from 10 to 100 MeV in increments of 
10 MeV, are the smallest values which will not be exceeded by actually 
encountered intensities. 

In order to hold the size and complexity of the code to a minimum, 
the desirable elements of the King (1974) analysis were empirically for- 
mulated by convenient approximation techniques, which are able to repro- 
duce the statistically treated King results with relatively high accuracy. 
Specifically, ti e maximum error introduced into the statistics by this 
method is less than 507.. This error is indeed insignificant in view of 
the large intrinsic statistical uncertainties of the models themselves 
and the much greater uncertainties involved in predicting cycle 21 fluxes 
from cycle 20 data. 

The most important feature of the code is the distinction between 
OR and AL events. This capability was developed relying entirely on 
King’s probabilistic treatment and using the data given therein ( King , 1974 
Figure 5). For simplicity let: 

n « f(Q, T) (I) 
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where n is the number of AL events expected to occur with a confidence 
level Q and for a mission duration t • It is then convenient to effect 
an apparent separation of variables by evaluating n independently for 
every value of Q within the regime of interest, and by expressing n 
in terms of T alone. Thus: 


6 

c + E 
j**0 


•ijC” 


( 2 ) 


where the index i ' relates n to Q through: 


i - 100 - Q (3) 

and where the polynomial expansion coefficients a.^ are determined by 
curve fitting. The constant, c=1.001, is a dimensionless range adjust- 
ment term. 


Table 1 contains the numerical values of the a 


ij 


coefficients. 


The code calculates n at the very beginning from equation 2 (see 
Table 2 for obtained values of n) and, depending on whether AL events 
are expected to occur or not, it branches to Its appropriate section, 
where it computes the corresponding AL or OR solar proton fluences 
for the specified Q and T . 


For AL-event conditions, that is n>0, the fluences are obtained 
from the expression: 

f al ( >E! Q,t > *' S al ( >E) * n( Q»T > (4) 

where the function: 


S ( >E) - 

AL 


7.9 x 10 9 exp {(30.-E)/26.5 } 


(5) 
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Is King * s (1974) exponential-in-energy representation of the integral 
flux and spectrum of the August, 1972, anomalously large event, for 
energy thresholds between 10 and 200 MeV. 

In case OR events only are predicted (n=0), the respective fluences 
are obtained from: 


f or ( >E;Q, t) - S QR ( >E)*J(Q, T) (6) 

where now King 's (1974 J expression: 

S. R ( >E) = exp (. 0158(30. -E) ) (7) 

for the integral energy spectrum of ordinary events is used, reliable 
over the range 10<E(MeV) <100. 

The function J(Q,t ) gives the ordinary fluence of protons above 
30 MeV in terms of the two input variables (see Table 3 for obtained 
values of J). It was developed along similar lines as equation 1, with data 
from the sane source ( King , 1974; Figure 8). 

Again, through the independent evaluation of J, (but this time for 
every r within the relevant Interval of 1-9 months for these events) the 
same apparent separation of variables is effected, resulting in: 

4 

J k “ d Z b kin (T)Q ra (8) 

m*0 

where the index m relates J to T through: 

m ■ integer(T) (9) 

and where the coefficients, contained in Table 4, are also determined 
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by curve fitting. The constant, 


d = 10^, is a scaling factor. 


The code is structurally complete and self-contained. All expansion 
coefficients for the statistical treatment are stored in data statements 
at the beginning of the routine, while the expressions for the AL and 
OR models are included in the respective sections of the code. 

To guard against accidental misuse and so as to insure that calcula- 
tions are only performed within the valid ranges of the independent vari- 
ables Q and T , S0LPR0 tests the values of these input parameters in 
its first executable statement. If the parameter(s) exceeds the program 
limits, it prints a warning message before returning control to the 
calling routine. 

Finally, the code bypasses ordinary event fluence calculations, 
although it branches to the appropriate section, whenever and Q 
has values ranging from 80 to 83. In tha - ' domain the flux gradients are 
very high, the intensities decrease very rapidly towards Q = 80 and the 
galactic cosmic ray component becomes predominant (see discussion by King , 
1974). For all practical purposes the OR-event solar protons are negligible 
for Q < 80 and t> !• 
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Resul ts and Discussion 

Sample calculations performed with the code indicate good agreement 
with probabilistic predictions. That is, the results represent a reasonable 
approximation to these data, considering the large intrinsic limitations 
and uncertainties inherent in this type of estimation process. 

Thus, as shown in Figure 1, the ordinary event fluence J of protons 
with energies E > 30 MeV, computed by the routine for the applicable Q and 
T ranges, compares well with the original data in King (1974; Figure 8); 
the solid curves indicate the former, the dotted curves the latter. Appa- 
rently, the curves overlap for most Q and t values. In the few places 
where they do not agree, the differences are truly insignificant; for 
example, the average error for all t> 1 is about 37., while the maximum error, 
which occurs at T=1 for Q=85 Is less than 507.. 

Similarly, the number of anomalously large events n given by the 
code agrees favorably with the model ( King , 1974; Figure 5). This is 
shown in Figure 2, where the dots indicate model data and the triangles 
S0LPR0 predictions. Since the transition from one n value to another is 
probabilistically (and not deterministically) defined, the differences 
between the King model and the S0LPR0 approximation are inconsequential. 

Finally, the AL event exponential-in-energy function (Equation 5) and 
the ordinary-event integral-energy spectrum (Equation 7) are identical to 
the formulae given by King (1974). Their respective values are listed in 
Table 5. 
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App 7 (cation and Use 


Subroutine SOLPRO is written In standard FORTRAN-IV computer language 
and card decks are available in either the 029 model IBM keypunch format 
(EBCDIC) or the 026 keypunch format (BCD), for use with FORTRAN-compatible 
compi lers. 

The variables in the calling sequence are: 

SOLPRO(TAU, 10, F) 

The first two arguments are input data: 

TAU : mission duration in months, with a valid range 
from 1. to 72. (REAL*4) 

IQ : percent confidence level that calculated fluence 
will not be exceeded, with a valid range from 
80 to 99 ( INTEC ER*4) 

The last argument Is output data: 

F(N) : solar proton spectrum for energies E >10*N MeV, 
where index N goes from 1 to 10. 

All aliments except IQ are single precision floating point variables. 
Besides their range limits, there are no other restrictions or limitations 
on the values of the input arguments. 

There are no READ or other input statements in SOLPRO. A WRITE 
statement is executed each time an input variable exceeds the stipulated 
limits. In addition; a WRITE statement is executed whenever for a given 
combination of T and Q no significant solar proton fluences are to be 
expected. 

9 




Figure 3 shows a flow diagram of the program. A complete listing 
is given in the attachment. The cards of the deck are appropriately 
labeled in columns 73-80 as SOLPRxxx where the three last columns (xxx) 
contain the sequential numbering, which is incremented by 10. 
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Table 1 
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NUMBE R Of ANOMALOUSLY LARGE EVENTS 
Polynomial Expansion Coefficients 


Q <%) 

•o 

*1 

*2 

*3 

*4 

*5 

*6 

99 

-.1571 

.2707 

-. 1269 “* 

, 4428’ 3 

-. 8185’ 5 

. 7754" 7 

-. 2939* 9 

98 

0 

CO 

• 

1 

.1951 

-. 6559 ’* 

. 1990’ 3 

-. 3618’ 5 

. 3740' 7 

-9 

-.1599 * 

97 

-.2007 

.1497 

-. 3179 ’* 

. 5730’ 4 

-. 4664’ 6 

. 1764 ’® 

0 

96 

-.1882 

.1228 

-. 1936 ’* 

. 2660 -4 

-. 1022 -6 

0 

0 

95 

-.2214 

.1149 

-. 1871 ’* 

. 2695" 4 

-. 1116’ 6 

0 

0 

94 

-.2470 

.1062 

-. 1658 ’* 

, 2367’ 4 

-. 9465’ 7 

0 

0 

93 

-.2509 

. 8710' 1 

-. 8300’ 3 

. 8438” 5 

0 

0 

0 

92 

-.2923 

. 8932* 1 

-. 1023 ’* 

. 1029” 4 

0 

0 

0 

91 

-.3222 

. 8648 -1 

-. 9992" 3 

, 9935" 5 

0 

0 

0 

90 

-.3518 

. 8417" 1 

-. 1000 ’* 

. 9956' 5 

0 

0 

0 

69 

-.3698 

. 7951" 1 

-. 8983‘ 3 

. 8940’ 5 

0 

0 

0 

88 

-.2771 

. 5473” 1 

-. 1543' 4 

0 

0 

0 

0 

87 

-.2818 

. 5072 ”* 

. 2511’ 4 

0 

0 

0 

0 

86 

-.2845 

. 4717 ’* 

. 5664’ 4 

0 

0 

0 

0 

85 

-.2947 

. 4405 ”* 

. 8507’ 4 

0 

0 

0 

0 

84 

-.2923 

. 4111 ’* 

. 1106’ 3 

0 

0 

0 

0 

83 

-.2981 

. 3853 ’* 

. 1312’ 3 

0 

0 

0 

0 

82 

-.3002 

• 3585' 1 

. 1529’ 3 

0 

0 

0 

0 

81 

-.3001 

.3312 1 

• 1781 -3 

0 

0 

0 

0 

80 

-.3141 

. 3248" 1 

. 1654’ 3 

0 

0 

0 

0 
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Table 2 

Number of Anomalously Large Events Predicted by SOLPRO 

99 
95 

| so 

85 
80 

*NO CALCULATIONS PERFORMED 


| I I | | I I I I I | I “ I 

— 1111122222222333333333334444444444555555555556606666667777777778 88 B8 8 888 
111111122222222223333333333333444444444444555555555556666666666777777777 
1)111111112 222222222223 33333333333334 44444444 4 4 44 555555555556666666 667777 
011111111112222222222222333333333333333444444444444455555555555566666666 

-011111111111222222 222222 22233333333333333344 444444444444 5555555555566666 - 
001111111111112222222222222223333333333333333444444444444445555555555566 
001111111111111122222222222222223333333333333333444444444444445555555555 
000111111111111112222222222222222233333333333333333444444444444455555555 
000111111111111111222222222222222222333333333333333333444444444444455555 

- 0 000111111111111 1112 2222 22222 2222222 22 33333333-35333333334444444444444555 - 
000011111111111111111222222222222222222223333333333333333344444444441.445 
000001111111111111111112222222222222222222333333333333333333444444444444 
000001111111111111111111222222222222222222223333333333333333334444444444 
000001111111111111111111112222222222222222222333333333333333333344444444 

-000000111111111111111111111222222222222222222223333333333333333333444444- 

000000111111111111111111111117.222222222222222222233333333333333333344444 

*00000011111111111111111111111222222222222222222223333333333333333333444 

*00000001111111111111111111111122222222222222222222233533335333333333344 

*00000001111111111111111111111111222222222222222222223333333333333333334 

-*00000000111111111111111111111111122222222222222222222233333333333333333- 

1111 I, 1 I I I I \ I I I 

10 20 30 40 50 60 70 

T (months) 






Table 3 


OR-Event Fluences (E > 30 KeV) given by SOLPRO Routine 


99 

98 

97 

96 

95 

94 

93 

92 

91 

90 

89 

88 

87 

86 

85 

84 

83 

82 

81 

80 


2-1 2 3 4 5 6 7 8 


5.07x10® 

3.42x10® 

2.34x10® 

1.67x10® 

1.30x10® 

1.09x10® 

9.78xl0 7 

8.85xl0 7 

7.72xl0 7 

6.20xl0 7 

4.34xl0 7 

2.40xl0 7 

8.86xl0 6 

5.23x10® 

* 

* 


* 


* 


4.83x10® 

3.91x10® 

3.23x10® 

2.74x10® 

2.39x10® 

2.12x10® 

1.91x10® 

1.73x10® 

1.56x10® 

1.38x10® 

1.19x10® 

1.00x10® 

8,16xl0 7 

6.56xl0 ? 

5.48xl0 7 


6.46x10® 

5.09x10® 

4.15x10® 

3.53x10® 

3.14x10® 

2.88x10® 

2.70x10® 

2.53x10® 

2.36x10® 

2.17x10® 

1.95x10® 

1.74x10® 

1.56x10® 


7.82x10® 

6.02x10® 

4.84x10® 

4.12x10® 

3.63x10® 

3.42x10® 

3.23x10® 

3.04x10® 

2.83x10® 

2.53x10® 

8 

2.32x10 


9.16x10® 

6.76x10® 

5.44x10® 

4.74x10® 

8 

4.32x10 

3.98x10® 

3.60xlC® 

3.22x10® 

2.98x10® 


8 

7.52x10 

6.25x10® 

5.39x10® 

4.80x10® 

4.38x10® 

8 

4.04x10 


8 

8.06x10 

8 

6.94x10 

6.08x10® 

5.41x10® 


l.OlxlO 9 

8 

8.38x10 

7.34x10® 


• Ho Calculations performed! see comments page 7. 
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Table 4 


ORDI NARY EVENT FLllENCE 
■Polynomial Expansion Coefficients 


(mths) 

b 0 

b l 

1 

. 154047xl0 3 

-.52225Sxl0 4 

2 

.198004xl0 3 

-.448788xl0 4 

3 

.529120xl0 3 

-.122227xl0 5 

4 

.121i41xl0 4 

-.266412xl0 5 

5 

.452062xl0 4 

-.103248xl0 6 

6 

. 272028xl0 4 

-.499088xl0 5 

7 

,275595xl0 4 

-.469713xl0 5 

8 

,570997xl0 4 

-.7996S9xl0 5 

9 

.101000xl0 3 

0 


b 2 b 3 b 4 


.714275xl0 5 

-.432747xl0 6 

.955315x10* 

.438148xl0 5 

-.196046xl0 6 

.325520x10* 

.112869xl0 6 

-.465084xl0 6 

.710572x10* 

.226778xl0 6 

-.857230xl0 6 

.120444x10 

.896085xl0 6 

-.346028xl0 7 

.499852x10 

.353050xl0 6 

111929x.l0 7 

.133386x10 

. 314729xl0 6 

-,960383xl0 6 

.111650x10 

.381074xl0 6 

-.610714xl0 6 

0 

0 

0 

0 
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Table 5 


Integral Flux & Spectrum 


of the August 1972 AL Event 


Distribution of 


OR-Event Integral Spectrum 


(Equation 7) 


Energy Solar Protons ( >E) 
(MeV) (#/cm^/AL-event) 


Spectral Function 
(normalized to 30 MeV) 


10 

1.680 

X 

10 10 

10 

1.37163 

20 

1.152 

X 

10 10 

20 

1.17117 

30 

7.900 

X 

10 9 

30 

1.00000 

AO 

5.A17 

X 

10 9 

40 

0.85385 

50 

3.71A 

X 

10 9 

50 

0.72906 

60 

2.547 

X 

10 9 

60 

0.62251 

70 

1.7A6 

X 

10 9 

70 

0.53153 

80 

1.197 

X 

10 9 

80 

0.45384 

90 

8.210 

X 

o 

00 

90 

0.38752 

100 

5.629 

X 

10 8 

100 

0.33088 




* amsaOM.* 
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Sample Program Listing 


SUBROUTINE SOLPP.O(TAtl, I Q,F> SDLPR010 
I NTERPLANETARY SOLAR PROTO!! F LUX AT 1 All (EROII £>10 TO E>100 MEV) SOLPR020 
SINOLE PRECISION BECK III STANDARD FORTRAN IV FOR IBM 360 MACHINES SOLPR030 
(EBCDIC, 029 PUNCH) OR OTHER COMPATIBLE SYSTEMS. SOLPR040 
PROGRAM DESIGNED AND TESTED BY E.G. 5TASS I I'OPOULDS, CODE 601, SOLPR05D 
NASA GODDARD SPACE FLIGHT CENTER, OPFENBELT, MAR, LAND 20771 . SOLPROEO 


INPUT: TAU 
1 0 

>•* OUTPUT: F ( N) 


MISSION DURATION IN MONTHS (REAL*4) 

CONFIDENCE LEVEL THAT CALCULATED FLUENCE F(N) 
WILL NOT nE EXCEEDED (INTFGER*4) 

SPECTRUM OF INTEGRAL SOLAR PROTON FLIIFNCE FOR 
ENERGIES E>10*N ( 1 -<!■-< 10) 


S9LPR080 

SOLPR096 

SOLPRIOO 

SOLPR110 

SOLPR120 


SOLPR290 
SOLPR300 
SOLPR310 
SOLPR320 
50LPR330 
SOL PP.34 0 


REAL HALE, NALECF ( 7, 2 0 ) /- . 15 71, . 2 70 7, - . 1269E-1, . 4 42 3 £-3, - . 81 8 5E-5, SOLPR130 
$. 7 754E-7, -. 2939E-9, -.1870, . 105 1, - . 6550 E-2, . 199 0E-3, 3618E-5, SOLPR140 

S.3 74CE-7,-.1500E-0,-.200 7, . 140 7, - . 31 70E-2, . 5 730 E-4, - . 4 664 E-F, SOLPR15 0 

$. 17C4F.-8, 0., 1802, . 1228, 1036E-2, . 266 0E-4, 1022C-E, 2*0 . , SPLPR160 

*-.2214, . 1149, -.1S71E-2, .2605E-4, -. 1116 6-6,2*0., -.24 70, . 1062, SOLPR170 

$-.1658 £-2, .23£7E-4,-.0465E-7,2*0.,-.2509, . 8710E-1,-. 8300E-3, S0LPR180 

*. 34 38 6-5, 3*0., - .292 3, . 393 2E-1, - . 1023E-2, . 1 029 F.-4, 3*0. ,-. 3222. SOLPRIOO 

$. 8648 E-l,-. 9902 E-3 , . 093SE-5, 3*0., -. 3518, .841 7E-1, - .1000F-2, SOLPR200 

$. 9956 £-5,3*0., -. 3608,. 795 IE-1,-. 808 3E- 3, .80406-5,3*0., -.2 771, SOLPR210 

$. 54 7 3 E-l,-. 15636-4,4.0., -. 281 8, . 50 72 E-l, . 2 511 E-4, 4*0 .,-. 284 5, SOLPR220 

*.4 71 7 E-l,. 56646-4,4*0., -.2947, .44056-1, .850 7 E-4, 4*0., -.2923, SOLPR230 

*.4 1116-1,. 1106 E-3, 4*0 . ,-.2081,. 38 53 E-l,. 1312 E-3, 4*0. ,-.30 02, SOLPR2 40 

*. 358 5E-1, .15206-3,4*0. , -. 3001, . 33 12 E-l, . 17816-3,4*0.,-. 3141, SOLPR25D 

$. 32486-1,. 16 54 E-3, 4*0./, F( 10), G( 10) SOLPP260 

REAL ORFLXC( 5,9 )/. 154 04 7 E3, • . 52225864, . 7 142 75 E5, -. 4 32 74 766, . 9553 15SOL PR2 70 
*E6, . 108004 63, - .4 48 78 8 C4, . 4 30148 65,-. 106046 66, . 32 5 52 E6, . 5291 2 0F3, BOLPP.2 8 0 
$ -.122227 65,. 112060 EG, -.4F 50 84 66,. 710572E6, .121141 64, -.26 6 412 E5, SDLPR200 

*.22677866, -.85723E6, . 12P444E7, .4 52062 E4,-. 10324866, .89608566, SOLPRIOO 
$-. 34602 8 67, . 4 908 5267,. 2 7202864, -. 4 0908 8 E5,. 35 30566, -.11192967, SOLPR310 

$. 1 33386 E7, . 27559764, -.4 F9713E5, . 314 720 E6, -.960383F6, .1116567, SOLPR320 

$.57009764, -.79068065. .30 10 74 £6, -.6 10 714 66,0., .10163,4*0./ 3DLPR330 

INTEGER INDEX(20)/2*7,6, 3*5, 5*4, 9*3/ SOLPP.340 

1 FORMAT ( ' TAN-*,F4.0, 1 I q- *, I 3, 3X, 1 PARAMETER(S) EXCEED PROGRAM lIM I 50LPR350 

$TS*> SOL PR 3 6 0 

2 TORMAT (2 X, 1 FOR THE COMBINATION OF TAU Ann 10 GIVEN, NO SI Gill F I CANTSOLPP.370 
$ SOLAR PROTON FLUXES ARE TO BE EXPECTED. TAU- 1 , F6 . 2, ' |D.',|2) SOL PR380 

lr(TAII.GT.72..OR.|O.LT.30)GO TO 500 SOLPR390 

I P-100- I q SOLPR400 

M-INDEX(IP) SOLPR410 

NAL6-0 . 5OLPR420 

DO 300 J-1,M SOLPR430 

0 NALE-MALC* NALECF ( .1, I P) *TAII**( J-l ) SOLPR440 

I NAL6-NAL6* 1.0001 S0LPR450 

IF( INALE. GT.O) GO TO 400 SOIPR460 

• CALCULATIONS FOR OR-EVENT CONDITIONS SOI.PR470 

I T -TAU SOLPR4SO 

IF( IT. 60.1. AND. IP.GT. 16) GO TO 700 SOLPR490 

P-F LOAT ( I P)/ 100 , SD1.PR5U0 

OF -0 . ‘ 5OLPR510 

DO 100 .1-1,5 SOLPR520 

0 OF -OF ♦ORFLXCL.I, IT)* P**( J-l)«l. E7 SOLPR530 

E-10. SOI.PRS40 

DO 200 11-1,10 SOLPR550 

G(N)-EXP( .0158*( 30.-E)) 50LPP.560 

F( N) -OF *G( II) SOL PR 5 70 

D 6-6*10 . SOLPR580 

GO TO 800 SOLPR590 

* CALCULATIONS FOP. AL-EVENT CONDITIONS SOLPR600 

1 E-10. SOLPR610 

DO COO 11-1,10 SOLPR620 

F( II) -7.9 E9*EXP( ( 30. - E)/26. 5)*l HALE SOLPR630 

0 E-E+10 . SOLPR640 

GO TO COO SOLPR650 

1 WR I TEC C, 2 ) TAU, IQ SOLPR660 

GO TO 800 SOLPR670 

9 WRITE (6,1) TAU, iq S0LPR680 

3 RETURN SOLPP690 

END SOLPR700 
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Sample Calculations: Sola* Proton Fluences F(>E;Q, t ) 



\m\ 


2 


7 


4 


r, 


6 


7 


n 


a 


10 


j 1 


1 2 

10 5. 7 l HC 

07 

2. r*6f 

44 

•». ’8*u 

0- 

1. 0* 

0 <3 

s .€ t ’t 

r i 

1 .AflCfc 

1 3 

! .6B4C 

l 0 

1 .c*HCE 

10 

1 .Cb3*L 

1 3 

1 #6ft3C 

10 

1 .68 OE 

l 0 

1 *6b3E 

X 0 

20 * • #w* ?t 

or 

1 » 79HF 

4 H 

1. 22 9t 

3 H 

% . 62 J« 

Ct 

h. .22 5C 

r h 

1 • 1 0?F. 

l 0 

1 • I bar 

1 0 

1 . 1 bar. 

1 0 

1 • 1 t>2f 

14 

1 • t 5? £ 

1 J 

1 .1 6 2fc 

1 0 

1 .1 S2h 

-X 3. - 

30 *• I6 1 * 

o r 

1 • 5 J4 F 

0 H 

2 • 75 7c 

3-. 

*. 1 lcJll 

: e 

7.32 Ifc 

' H 

r.«*Cfc 

CS 

7.9<*4fc 

30 

/.Q4 0t 

39 

7 • 9 0 0 F. 

09 

7,9007: 

09 

7.90QC 

OS 

7 . 90 JE 

09 

*0 J. 56 0f 

o r 

1 • 01 IE 

4 

2. 054C 

3h 

3. 01 61. 

o t 

? .SS6t 

OH 

5.4 i ?r 

Ci 

j.4 1 Tf 

0 9 

5.* 1 7 t 

09 

5 • * 1 7f 

0 9 

b.* 1 7C 

09 

b . * 1 76. 

OS 

b.* 1 ?t 

09 

50 J • 0 3 91. 

or 

1 • 1 2 r F 

OH 

?. 4»6L 

3e 

J. 0 3 3f. 

c« 

t . 1 2 It 

f H 

1.71 4F. 

09 

i.7»*r 

09 

0 .71 4fc 

09 

3.7l*f 

09 

3 . 7 1 * F. 

09 

3 .7 ME 

0 s 

1 . 71*E 

09 

60 2.59 ?t 

•7 7 

«, 6SSF 

4 T 

1 . 7|oE 

)3 

2. 5»j*t 

OF 

* .372£ 

CH 

H.*i* or 

OS 

2* b* 7t 

09 

? • b* 7fc 

09 

2 • 54 7i* 

09 

2.4*75 

09 

2.b* 7E 

0 s 

2 . 64 re 

09 

70 2. ? 1 6t 

o r 

H « Ui?t 

6 7 

l • *661 

C ;* 

1*91 

: f 

: .73 

( H 

1 • 74tE 

09 

! . 7 4 tr 

09 

1 . 74 6C 

09 

1 . 746t 

4U 

I • 7*6 1 

09 

1 . 7*bt 

0 9 

1 • 7*t»E 

09 

8Q 1 • 89 2L 

3 7 

t • sr>9-‘ 

4 7 

1 . 26 1 t 

Or* 

1 . »6 9t 

n * 

1 . 1 ‘i rr 

C«5 

I . I SOE 

Os 

1 . IS7F 

09 

1 . 197fe 

09 

1 • I9rf, 

40 

1 • 1 S?l 

09 

1.1975 

0 9 

1 • 1 37£ 

09 

90 1 • 6 l V 

37 

h • *>*.> 1 t 

4 7 

! . 26‘JL 

M 

J • V»*>t 

Of 

2.72 !c 

* H 

•1 • :m Ct 

CH 

a.2i<3t 

3 a 

*»• ?1 3E 

4(3 

ft .21-37 

08 

b .2t 4 1 

Oft 

8.2 105 

3 8 

8.21 JE 

33 

100 1 • 17 9t 

o r 

b.^btE 

47 

>. 1? U 

■ ) 7 

1 . 

Of 

£ . J2 4t 

rn 

0 * 0 2 s t 

ca 

b. 6 29F 

3J 

4.629C 

Oft 

t»,629r 

3ft 

6.6 2Cc 

Ob 

5.629L 

00 

5.629E 

OH- 

>E(Ke?) 1 " 7 T*I J 


14 


!'• 


1 4 


1 7 


1 8 


to 


20 


21 


22 


2 3 


24 

10 i *o« *>. 

1 0 

1 »*tt6£ 

1 4 

1 .6H3* 

1 ' 

1 • 6 9 Ofc. 

1 * 

1 .c « je 

l 0 

1 .48'T 

1 0 

1 .4H6f 

t 4 

1.680C 

10 

\ .684C 

10 

1 .6805 

to 

1 .6ft0£ 

1 0 

1 .680b 

1 0 

20 l • 1 b 2h 

1 6 

1 • t*»?F 

1 4 

1 • l*2f 

1 

1 . 1 62l 

1 3 

1 . l'j.’r. 

l ' 

1 . lH?t* 

1 3 

1 • 1 52f 

1C 

1 . lb2E 

1 0 

1 . i52r 

14 

1 . 1 sac 

1 3 

1 .1 52£ 

1 0 

t . l 52C 

1 0 

30 7, 90 

4 9 

7, q^OF 

4 C 

7 • 900L 

0 * 

fm VOCE 

; < 

7 . SC 31 

c«- 

r. *:cc 

CS 

7.44 0t 

39 

7. 9 0 0 C 

39 

7.940T 

49 

7.903E 

>9 

7 .9 0 OE 

OS 

7. V03c 

09 

hO - • * 1 7b 

0 7 

••• 4 1 7F 

4 g 

6 . A | 7l 

4'. 

r ». *17. 

3<; 

t . 4 1 7P 

0 •» 

5.4 1 re 

OS 

S. 4 lT t 

09 

s.* 1 re 

09 

6 • * 1 7r 

OO 

?» . * 1 7 1 

39 

5.* 1 7c' 

0 s 

5 • * 1 7£ 

39 

50 • • 7 1 *C 

09 

J. 71 4f 

r q 

?. 7|*l. 

Ow 

1. 71 4r 


3 . 71 %l 

2 v 

3. ri4t 

OS 

J. 7 I4f 

09 

J. 71 4f. 

ng 

J. 71 4F 

49 

3. rut 

39 

3.7 1*5 

OS 

J. 7 J *6 

09 

60 <** 5* 7F 

09 

2 «4 ?,* 

<H 

2 • •?* 7r 

0 . 

2. *4 7. 


2 .t* 7*: 

C » 

? *64. 7F 

eg 

?.447f: 

49 

2.S47E 

09 

2.4* 7L 

09 

2.S* 7C 

09 

2 .6* 7c 

0 s 

2 . 6 4 7b 

0 9 

70 1 • 7* tx- 

09 

1 • 74 f £ 

4q 

1 . 746c 

O'* 

1 • , 4 4.. 

r s 

1 . 74* * 

' > 

l .74ft 

09 

1 .74f.r 

49 

1 . 746t 

49 

1 . 74 t>r 

09 

1,7 tot 

CO 

1 .7*65 

0 s 

1 • 7*oE 

0 9 

80 1 • l'77t 

09 

! • 1 9 7 r 

09 

1 • lorn 

3 r * 

1 . I97r 

4 <; 

1 . » S 71 

c>; 

1 • I S T* 

39 

1 . 1 97r 

09 

1 . 1 07F 

09 

1 • 1 9’r 

49 

1 • 1 S7E 

09 

1 .1 >te 

OS 

1 • 1 9.‘E 

09 

90 8. 2 1 0* 

Ort 

ft. Jl rt t 

4 rt 

H,2|4t. 

3y 

0*21 0 1. 

0 1 

* . 2 1 1L 

c 0 

n.i'i cfc 

OH 

H.2 1 3E 

OH 

H.21F t 

Oft 

8.21 0£ 

Oft 

4 .21 >e 

08 

8.21 )C 

0 e 

8.210b 

0 8 

100 5 . oJ 9, 

0*5 

t. 6? Of 

6p 

b • * 2 6 c. 

0“ 
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Of 

t . 7. ? 9t 

C H 

6.4?Sl 
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b.o2 

Orf 
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Ob 
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Ob 
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OF 
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1 1 

i. »6lf 

1 : 

2 • j 4 1 r 

1 > 

3 * Jo IK 

l - 

1. IMP 

1 0 

1. 30 IF 

1 0 

J . 1417 

i 0 

J.i6*r 

10 

J. 36 IE 

1 0 

3.361b 

1 0 

20 2 • 3.* «•«. 

1 6 

?• ^* r 

1 4 

2^ *44 L 

l 

2. «**t 

i : 

2 . 1 )1C 

I ' 

2 . 1 0 *C 

1 : 

?* 1 4* *3 

1 0 

2 ,ir*F. 

10 

2. 16*1 

1 0 

2.i o*r 

1 0 

2.3 3*c 

1 0 

2.3045 

1 *3 

30 i • *h oe 

l C 

1 • **f p 

1 4 

t . *501 

I* 

l . ooor 

1 ? 

t .bdOK 

1 *» 

1 .*>*jrr 

1 4 

1 .4?4 L - 

10 

1 • bH 3 F 

1C 

\ .seor 

10 

t .6H0c 

10 

l .OttOt 

1 c 

l . j»OC 

1 0 

40 X • Cb Jh 

10 

1 . or 

1 0 

l . CH Jk 

1 . 

1 . :•> :f 

1 3 

1 .CH »l 

1 3 

1 .■»!« ?r 

1 0 

1 .OF JL 

l 0 

1 .Ob JF 

10 

1 *0ft3 r 

10 

1 .Oblfc 

10 

1 .Obit 

1 0 

1 • 9b3£ 

1 0 

50 7, *2 Hfc 

09 

7, 4?bf 

4 9 

7. 4?HL 

' (' > 

7.4 .*«•. 

;s 

7 .42 5C 

C 9 

7 . * 2 f r 

•39 

7.42HC 

09 

7. *?HE 

0 9 

7.42UC 

09 

/.* »Hh 

99 

7.428E 

0 s 

7.A2BE 

09 

605,09 31 

0 9 

6. 09 jr 

4 9 

'•.49 »r 

4 9 

•>.C 13c 


3 . C < J F. 

: 1 

6 • "'S 2f 

0 9 

4. ?9lf 

4 9 

8.09 Jt* 

09 

b .09 ir 

09 

6.09 JL 

09 

6.99JE 

3 S 

S.0 9 i£ 

09 

70 J • * * 21 

J * 

:• 4 >2* 

4g 


',•/ 

J. * 72* 

39 

2 .*9jr. 

0 V 

J.4S ?f: 

C 9 

1.* V?t 

00 

1.492f 

09 

J.*9?F 

no 

J .*91C 

09 

3«*92C 

0 S 

J.«92C 

09 

80 ?• 3-vbe 

09 

2. 195t 

4 g 

?• W‘JL 

r. w 

2. Ubfc 

C9 

2 . 39 3t 

n 1 

?. 5S “C 

39 

2. 3V>C 

39 

2. J95E1 

0 9 

2. 396F 

09 

2. J9bfc* 

09 

2. JVot 

0 s 

2.3 V5£ 

09- 

90 1 • 64 a. 

39 

lit *2» 

49 

1 • ' *21 

3> 

1 , r.M 2L 

09 

t . 1 4 .;r 

C ) 

1,6 4 ■»F 

OS 

t .4*2F 

09 

1 .6*?E 

39 

1 .6421 

69 

1 .6*2t 

30 

1 .6* 2b 

OS 

1 .6*2£ 

09 

100 X • 1 2 St 

09 

1 • 12*1 

4 q 

1 • I26u 

Ov 

1 . I Jot 

os 

1 . 1 '6fc 

Z 9 

1.12th 

09 

» ♦ 1 ?6E 

OS 

l .1266 

09' 

l . 12bt 

09 

1 • 1 2o»: 

09 

1.126C 

OS 

1.1 26b 

09 
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4 0 

10 3. 36 lc 
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2. "*6 1 c 

1 4 
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i< 

M.k 

1 3 

2 . 14 ir 
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1. l6tt 

1 0 

1,.l6ir 

1 0 

J. Jolt 

10 

b.o*ir 

1 0 

6 • 0* 1 r 

10 

5.o*in 

1 C 

5.0*1 E 

1 0 

20 2m J3*i 

l C 

2. 14**. 

1 4 

.*04L 

1 ' 

.■*. 1 '4t 

1 r 

2 . 104t 

1 3 

a. 104F 

1 0 
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1 0 
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10 

3.*bot 

10 
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10 

3.*bnt 

1 0 
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30 i . t»c5 oi 
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1 1 

l.ifk 

J4 
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1 r 
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10 
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1 .62 4F. 
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09 
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09 

7 • *2bt 

09 

1 . u*r 

14 

1 .1 1*1 

’.0 
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*: • c 9 if. 

C > 

6 . 3 S It* 


6.0 Jlh 

40 

4.49 IE 
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09 

70 J. **2t 

09 
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-q 
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09 
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09 
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b.2 JVE 

09 
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l.6-*2c 

39 

1 . C 4 2t 

C 9 

t • 6* ?r 

09 
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09 
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09 
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09 
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09 

1 .7>B Vfc 
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09 
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09 
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H.4* 1 C 

1 4 

6. Mk 

1 ' 
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10 
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10 
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3.46tC 
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1 c 

3.456E 

l 0 

30 2* 3 7 CL 

1 ( 

2 , 17^f 

l 4 

2. ’70L 

1 • 

2 • 1 T Oi 

1 : 
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10 

2. J70r 

14 

2.370c 

1 0 

2>iTdE 

1 c 
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1 .1 1*F 

’l 4 

1.1 1«( 

1 0 

1.11*1 

l 0 

1.1 l*h 

14 

1 • t 1 * c 

1 3 

1 .1 t*E 

1 0 

l • ' 14b 

1 0 

60 7.64 Ol 

OV 


r g 

7 . t>* 3 f 

4^ 

7. f»4 0t 

o« 

7 .64 3C 

C 7 

7.04 (r 

OS 

7 . t * 3t 

09 

7. t»*or 

09 

7.640L 

09 

7,0* or 

39 

7.6* Of. 

09 

7- ,4.0S 

OH- 

70 b.219: 

09 

6.2 *9‘_ 

r 9 

2 IV fc . 

s'V 

0 • 2 Wt 

c s 

? .23 ;c 

C V 

4 . 0 JST 

49 

6.2 ivr 

09 

b, 2 JSC 

09 

b .2J9t 

49 

S ,2 19 c 

0 t 

5.2 19C 

OS 

6.239E 

39 

80 '.b9?*' 

■>« 


r 9 

4.b92c 

3 1 

6 7?t 

3 S 

2 » 2c 

4 1} 

1 • b 0 2 E 

OS 

». e ;’c 

09 

1 . t> 92 t 

09 

J.S92* 

09 

j.boar 

09 

3.*» V2L 

OS 

J.592E 

39 

90 2. 40 3t 

0 u 

2.46 ir 

4g 

2* 46 JL 

06 

2. *4.17 

OS 

2 .4t 3* 

C 9 

?.4f» Jfc 

**s 

.* • * * IE 

09 

2 • *t>2 r 

09 

2 . * 6 1E- 

4q 

2.461F 

09 

2.463E 

OS 

2.*63E 

09 

100 1.6***- 

3v> 

1 .6ft9t 

4 Vi 

1 • t'BvE 


t .6-3 9c 

49 

1 .C3S»: 

f -> 

» • 4H SC 

39 

1 . •> H )F 

00 

1 .689C 

09 

1 .6isr 

09 

I .6b9t 

39 

1 .639E 

OS 

X -699E 

09 

>E(He?7 f trel 


* 2 




C 4 


f 


64 


S 7 


68 


64 


70 


7 1 


72 

10 *. '* 

J V • 

!»• '*n 

1 4 

6.72U 

iV 

6. 7*.l c 

1 C 

6 , 72 ir 

1 3 

6.721C 

10 

7. • 7 2 1 f 

1 0 

6 . 72 U 

10 

0 * 72 1 r 

14 

6 .721 r 

10 

6.7aifc 

1 0 

6.721 E- 

13- 

20 1.**j6c 


2* 4* 6 r 

1 4 

4 . 

1 * 

4.600, 

1 c 

4 ,o3 9c 

1 ** 

4 .4 osr 

l 3 

* . 4 J »t 

\ 4 

4 .t>09F 

1C 

* ,60 9? 

10 

* ,6-» >r 

10 

* .03 9f. 

l 0 

* .609b 

10 

30 2. J7Cc 


2. 47*r 

1 * 

'♦ 16 3t 

14 

1. loot 

1 c 

2 • 16 3C 

I ' 

1.1 6CE 

1 3 

j.i 4^r 

1 0 

1. lft f t 

to 

■». 16 OK 

1 2 

1 . !• 't 

10 

3.1 60 t 

1 c 

3 . 1 60E 

1 0 

40 1 • 6 ? V 

l 0 

1 . *2«-r 

1 4 

14 71 

1 : 

?. 14 7r 

1 c 

Z . 1 4 7c 

! 0 

2.14 ?t 

1 3 

' . 1 r> 7t- 

1 0 

2. I0 7C 

10 

2.17' 7 r 

14 

2. 16 71 

1 3 

2 . 1 6 7F. 

1 0 

2. 167E 

1 0 

50 1*1 14*- 

1 J 

1. U4r 

14 

1 . 4f*i»C 

r 

1 . 4 ,4c 

1 0 

1 • 4H6t 

1 ** 

1 • * r t h 

l 3 

1 .4 a *>r 

1 0 

1 .497 L 

10 

1 .4 84F 

14 

l .* 14 t 

10 

t ,*b6 £ 

l c 

1 . * H6b 

l 0 

60 7 • o* 3 l 

09 

7. ** ft 

4 Q 

1 . 4 1 >t 

1 ‘ 

1 . •*! 9t 

1 “ 

1 .’’l 7*. 

1 3 

l . ' ISf. 

1 4 

1 . 4 1 »L 

10 

t .0 19£ 

10 

1 .6 | If 

10 

1 .4 Wfc 

1C 

1.319c 

1 0 

1 • 3 1 v£ 

1 3 

70i.2JS* 

4v, 

b. 2 1st 

4 9 

r .C4t,c 

j-* 

cv. '.Hbc 

0 s 

c . <. t oc 

cs 

6. if*bf 

C J 

4. 

09 

6.98f»L 

4«J 

6.956C 

49 

6. 

09 

c.gftbt 

OS 

0 . 9*J6b 

09 

80 J.5 <A 

39 

„*• ‘♦9?r- 

4q 

<■ • 73 n 

? / 

4. 7-»‘.r 

OS 

4 • f'i St 

C i 

4 . 7 e* sp 

39 

4 • 740f 

39 

* • 7P0c 

09 

4.7ft 3F. 

09 

* • TrtOf, 

>9 

* • 7 ft 9c. 

OS 

4 , 78 JE 

09 

90 V • 4*» If 

O'. 

• 4f. If 

4 9 

*• i'P*C 

'c 

». ?3*r: 

c s 

: .2*4f 

C 9 

1, .*«*r 

OS 

i. 2 «*r 

09 

J..*H*t 

0 9 

J . 2847 

49 

i . ? H* r 

39 

3.2 ft*t 

OS 

3. 2 b* 6 

09 

100 l . O0 9L 

09 

! . 6ft 9t. 

49 


fl, 

2. 2 V 21 

OS 

* . ? ^ 3L 

( 9 

2.252: 

C-i 

2.2* ;L 

C9 

2 ,252f 

39 

2 • J>?r 

49 

2 .2521* 

29 

< .252c 

OS 

2.2525 

09 





